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ERROR  STATISTICS  FOR  ASTROGEODETIC 
POSITIONS  FOR  AN  RGSS  TEST  COURSE 


INTRODUCTION 


The  objective  of  this  research  is  to  provide  error  statistics  for  the  series  of 
astrogeodetic  positions  in  a  test  course  required  to  evaluate  the  U.S.  Army  Engineer 
Topographic  Laboratories  (ETL)  Rapid  Geodetic  Survey  System  (RGSS). 


The  National  Geodetic  Survey  (NGS)  of  National  Oceanic  and  Atmospheric 
Administration  (NOAA),  U.S.  Department  of  Commerce,  observed  and  computed 
the  astronomic  latitudes  and  longitudes  for  the  ETL  test  course.  The  latitudes  and 
longitudes  were  obtained  at  each  station  of  the  test  course  from  observations  taken 
on  two  nights  with  a  Wild  T-4  Universal  theodolite  and  Datametrics  model  S-P 
300  digital  timing  system. 


Latitude  determinations  were  made  using  a  modified  version  of  the  Stcrneck 
method.  Longitude  determinations  were  made  using  the  meridian  transit  method. 
Time  synchronization  for  the  Datametrics  time  system  was  maintained  from  radio 
signals  transmitted  by  the  National  Bureau  of  Standard  Time  Service  Station  (WWV ) 
located  at  Fort  Collins,  Colorado.  Stellar  positions  taken  from  the  Fourth  Funda¬ 
mental  Catalogue  (FK4)  were  used  for  computing  latitudes  and  longitudes.  For 
latitude  determination,  the  modified  Sterneck  method  involves  two  zenith  distance 
measurements  of  both  occular  positions  that  are  symmetrically  east  and  west  of  the 
meridian.  The  time  of  the  stellar  bisection  is  recorded.  The  latitude  determination 
depends  on  observations  of  32  stars  with  occurrences  divided  north  and  south  of  the 
zenith.  Stellar  positions  taken  from  the  Fourth  Fundamental  Catalogue  were  used 
for  reducing  astronomic  latitudes  and  longitudes. 
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ANALYSIS  OF  ERRORS  IN  LATITUDE  DETERMINATION 
BY  THE  STERNECK  METHOD 


ANALYSIS  OF  ERRORS  IN  LATITUDE  •  Erroneous  latitude  deter¬ 
minations  result  from  observing  zenith  star  angle  errors  that  are  due  to  personal  equa¬ 
tion  of  bisection,  index  error,  refraction  uncertainties,  and  scale  reading  errors  from 
inaccurate  graduation  line  coincidence. 


The  index  and  bisection  errors  are  cancelled  when  two  observations  are  made 
in  the  direct  and  reverse  theodolite  position.  In  the  modified  Sterneek  latitude  method, 
each  star  is  observed  symmetrically  east  and  west  of  the  meridian  plane.  Corrections 
are  applied  to  the  observed  zenith  distances  and  reduced  as  the  zenith  distances  are 
observed  in  the  meridian  plane.  Thus,  the  latitude  can  be  determined  from  the  general 
formula: 


<t>  =  5  ±  Z  (1) 

where  5  is  the  star’s  declination,  and  Z  is  the  mean  value  of  observed  zenith  distances 
after  applying  all  the  corrections.  Latitude  is  determined  from  observations  of  32 
stars  with  occurrences  distributed  north  and  south  of  the  zenith. 


Evaluation  of  Personal  Error.  After  reducing  the  star  zenith  distances 
from  direct  and  reverse  instrument  position  to  a  common  plane  of  references,  one  can 
compute  the  index  error.  Let  Z,  be  the  observed  zenith  distance  in  the  direct  position 
measured  0°  -*  to  90°  and  Z2  be  the  zenith  distance  on  the  reverse  instrument 
position  measured  0°  to  270°.  Let  e  be  the  instrumental  index  error  that  results 
from  observing  a  star.  The  personal  error  is  small  and  can  be  treated  as  an  index  error 
because  the  observer  sees  the  star  before  or  after  the  star  crosses  the  horizontal  reticle 
line.  The  total  effect  of  the  instrument  index  error  and  personal  error  is  considered 
as  (t  +  p). 


A  fixed  object  with  direct  and  reverse  zenith  distance  observations  gives 
the  following  equation  of  condition: 

Z,  +  Z2  +  2(e  +  p)  =  rr  (2) 
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where  Z,  refers  to  the  direct  observation  and  Z2  to  the  reverse  instrument  position. 
The  value  of  p  differs  when  the  observation  is  made  with  respect  to  a  moving  object, 
and  (e  +  p)  must  be  referred  to  as  star  observations.  Owing  to  the  star  motion 
within  the  zenith  distance  observations,  corrections  are  to  be  introduced  to  the  values 
Zj  and  Z2  to  satisfy  equation  (2).  This  can  be  made  by  reducing  Z,  arid  Z2 
to  a  common  instant  of  time,  T0.  Let  T,  and  T2  equal  the  references  times  that 
correspond  to  the  observation  of  and  Z2  respectively.  Applying  corrections 
A7-!  and  AZ2 ,  reduce  Z,  and  Z2  to  the  time  T0.  But  to  fit  equation  (2),  one 
must  include  the  level  and  refraction  correction  also.  Let  I*!  and  R2  equal  the 
refraction  corrections  to  Z,  and  Z2,  respectively.  Therefore,  the  level  corrections  are 
p/2  b|  and  p/2  b2.  After  adding  the  corrections,  equatio;  (2)  changes  to 

Z,  +  Z2  +  R,  -  R2  +  (bi  -ha)  +  A2i  +  :<e  +  P)  “  2tt  (3) 

P  =  Level  value  in  seconds  of  arc  per  division  of  the  vial 

Hast  star  provides  an  equation  where  AZt  and  AZ2  are  known  through  computation 
and  the  personal  equation  is  the  only  unknown  to  be  considered.  If  n  is  the  number 
of  observed  stars,  the  personal  equation  is  derived  from 

P  =  ff"2^I£(Z|  +Z2)  +  I(Ri  -R2)  +  4(~(bl  “b2>  +  (2AZi  -SAZ2)]  -e  (4) 

EVALUATION  OF  AZ(  and  AZ2  •  Let  Z0  equal  the  zenith  distance 
that  corresponds  to  a  time  T0,  and  let 

T  =  T  -  T 
1 1  1 1  1  o 


According  to  Taylor’s  theorem. 


From  the  condition 


and  from  equation  ((>)  is  obtained 


sin  Z  — —  =  cos  6  (12) 

dt 

in  which  Z  represents  the  zenith  distance  on  the  meridian  plane.  If  a  star  has  an 
azimuth  A  =  0,  then  dA/dt  has  a  positive  value.  If  the  azimuth  is  A  =  it,  dA/dt 
lias  a  negative  value. 


Analysis  of  the  observed  field  data  shows  first  that  no  index  error  was  evaluated 
and  second  that  the  star  zenith  distances  were  not  measured  with  respect  to  the  central 
intersection  of  the  vertical  and  horizontal  lines  of  the  reticle.  The  theodolite  position 
was  fixed  with  respect  to  the  meridian  plane.  The  stars  were  collimated  outside  of  the 
vertical  center  line.  Specifically,  the  analyses  shows  that  the  index  error  disappears 
when  averaging  the  direct  and  reverse  zenith  distances.  Equation  (5)  cannot  be 
applied  to  obtain  the  zenith  distances  corrections  because  the  equations  were  derived 
for  observations  made  for  the  intersection  of  the  vertical  and  horizontal  center  vertical 
lines.  To  apply  equation  (5),  a  meridian  projection,  PZ,  is  used  on  the  horizontal 
plane  (figure  1). 


FIGURE  1.  Observing  a  Northern  Star. 
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/here  SMS'  is  the  parallel  of  the  north  star’s  declination.  Let  r  equal  the 
hour  angle  of  the  star  when  the  star  is  collimated  at  S.  Since  the  theodolite  is  clamped 
with  respect  to  the  meridian,  the  horizontal  line  of  the  collimation  is  SN.  The  zenith 
distance  ZN  is  measured  instead  of  the  zenith  distance  ZS.  With  meridian  zenith 
distance  as  ZM,  the  correction  is 

MN  =  PM  -  PN  =  PS  -  PM  (13) 

The  value  of  PN  is  obtained  from  the  spherical  triangle  PNS  as  a  function  of 
the  hour  angle  r,  and  the  polar  distaiv  e  PS  =  90°  -  6.  Setting  PN  =  90°  -  6', 
we  obtain 


tan  6'  =  tan  5/eos  r 


S'  is  very  close  to  6  and  allowing 


cos  r 


-  2  sin2 


2 

2 


T 


and 


6Z  =  6 


-  5  = 


sin  26  sin2  r 
sin  1 


"therefore, 


(14) 


MN  =  6  ’  -  6  =  6Z„ 
The  meridian  zenith  distance  of  the  north  star  is 


Zm  =  ZM  -  6Zn  (i 5) 

Where  SZ  is  always  subtracted  from  the  observed  zenith  distance  of  a  north  star. 
The  correction  to  be  applied  to  a  southern  star  is  shown  in  figure  2.  where  ZN  is  the 
observed  zenith  distance  and  the  meridian  distance  is  ZM. 


Thus, 


ZM  =  ZN  +  NM 


For  a  south  star,  tiie  observed  zenith  distance  is  less  than  the  meridian  zenith  distance. 
A  correction  6Z  -  6  '  -  6  is  added  to  the  observed  zenith  distance. 


p 


FIGURE  2.  Observing  a  Southern  Star. 

OBSERVATION  ERROR  •  The  Sterneck  method  yields  the  latitude  from 
the  generalized  equation 

0  =  6  ±  Z  (16) 

in  which  Z  represents  the  meridian  zenith  distance.  The  positive  sign  is  for  the  south 
star  and  the  negative  sign  for  the  north  star. 

In  the  NGS  modified  Sterneck  method,  the  index  error  is  eliminated  and  the 
|  need  for  the  zenith  point  correction  is  avoided.  The  mean  zenith  distance  of  a  star 

is  corrected  for  refraction  and  inclination  errors,  and  is  reduced  to  the  meridian  to 
;  obtain  the  Z  value  that  appears  in  equation  (16).  The  probable  error  of  an  obser- 

[  vation  may  be  found  by  comparing  the  values  of  tlie  index  error  and  the  personal 

I  equation  that  results  from  the  direct  and  reverse  star’s  zenith  distances.  For  each  star, 

|  the  observed  values  are  reduced  to  a  common  instant  of  time  after  applying  ali  the 
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corrections.  According  to  the  theory  of  least  squares,  where  e  denotes  the  error  of 
an  observation,  V,  denotes  the  residuals  obtained  by  comparing  the  mean  (e  +  p) 
of  the  first  star  with  n1  observations 


•  t 
\ 


n 


n 


m 


then, 


(nj  -  1)  ee  *  p2  [V,  V,  1 
(n2  -  1)  ee  =  p2  [V2  V2  ] 
(n,„  -  Dee  =  p2  [Vm  VmJ 


m 

4 


4 

I 


i 

* 


I  where  (V,  VJ,  IV2V2),  etc.  denote  the  sum  of  the  squares. 


The  sum  of  this  equation  yields 

(n  -  m)  ee  =  p2  [V  VI 


Ja 

h.; 

ks 


if- 

I? 

Kv 

f 

u 

? 

#;* 


where  n  denotes  the  whole  number  of  an  individual  star’s  zenith  distance  and  re¬ 
sults  in 


Hence,  we  have 


n  -  n,  +  n2  + 


+  n.. 


IV  Vi  =  IV,  VJ  +  [V2  V2]  + . +  (Vm  Vml 


The  NCiS  modified  Sterneck  formulation  involves  two  symmetrical  eircummeridian 
zenith  distances  in  direct  and  reverse  instrument  positions: 


I 


anti 

i 

'  in  =  n 


Hence,  we  have 

e  =  =p 


PROCEDURE  FOR  OBTAINING  OBSERVATION  ERROR  •  As  de¬ 
fined  in  the  previous  section,  the  values  of  V  come  from  the  star’s  zenith  distances 
on  the  meridian  plane  at  the  instant  of  time  T0.  Each  star  is  observed  from  the 
meridian  at  an  hour  angle  r  defined  by 

r  =  0o  +  T0(l  +C)  -  (X+«)  (17) 

where 

T  0  -  Greenwich  Universal  Time 

0()  -  Greenwich  hour  angle  of  vernal  equinox 

C  =  Constant 

X  =  longitude 

a  =  star  right  ascension 

The  r  values  given  by  equation  (17)  are  used  to  obtain  5Zy.  Let  Z,  represent 
the  direct  observed  zenith  distance  at  the  instant  of  time  T j ;  Rj  the  error  due  to 
refraction:  6Z(  the  correction  to  reduce  the  observed  zenith  distance  to  the  meridian; 
(e  +  p)  the  index  and  personal  error  of  observation;  and  Z  the  true  value  of  the  star’s 
zenith  distance.  For  direct  instrument  position,  the  following  equation  must  be  satis¬ 
fied: 

Z  =  Zj  +  Rj  +  (e  +  p)  +  7  •  5Zj  +  Vj  (18) 

where  y  symbolizes  the  sign  to  be  applied  because  the  correction  5Zj  and  V, 
is  the  error  of  observation.  A  similar  equation  is  obtained  for  the  reverse  instrument 
position,  which  is  indicated  by  the  subindex  2. 

Z  =  2jt  -  Z2  -  R2  -  (e  +  p)  -  7  *5Z2  -  V2  (19) 


Let 


Zj  Zj  +  Rj  +  ~  bj  —  7  •  5Zj 


Zr  =  2?r  -  Z2  -  R2  -  7  •  6Z2  +  —  b2 


(20) 


introducing  the  values  in  (20)  and  subtracting  (19)  front  (18)  results  in 

0  =  Zd  -  Zr  +  2(e  +  p)  +  V2  +  V,  (21) 

Each  star  furnishes  an  equation,  such  as  (21),  with  constant  value  (e  +  p).  With  n 
number  of  stars  there  are  in  equations,  which  when  added  give 


S(Z„  -Zr)  +  2n(e  +  p)  +  S(V,  +  V2)  =  0 
Solving  for  the  unknown  (e  +  p)  yields 

S(Zd-Zr)  SfVj+Vj) 


e  +  p  =  - 


2n 


2n 


The  V’s  are  random,  and  when  n  is  large,  then 

s(yi+v2>- .  o 

2n 


and 


e  +  p  =  - 


:(zd  -zf) 

2n 


(22) 


(23) 


(24) 


Knowing  (e  +  p),  the  V’s  can  be  defined  from  equations  (18)  and  (19). 


As  an  example  of  computation,  we  took  the  observation  of  Station  Aero 
USAETL  1978,  shown  in  table  1.  Table  2  is  self  explanatory.  The  hour  angles 
were  computed  using  equation  (17).  The  5Z’s  from  equation  (15)  and  the  re¬ 
traction  correction  from  the  author’s  formula  results  in1 


1  A.  A.  Baldini,  "Formula*  for  Computing  Atmospheric  Refraction  for  Objects  Inside  or  Outside  the  Atmosphere.” 
GIMRADA  Research  Note  8.  1963.  AD -419  915. 
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TABLE  1.  Modified  Stcrneck  Latitude  Computation 

1  STATION  AERO  USAETL  1978  MD  INST  NO.  T-4,  86994  CHRON  NO.  D-253 

2  OBS  BY  COHEN  REC  BY  JCH  ACC  NO.  ECC  O.OOOM  AZI  ODEG  O.OMIN 

3  LATITUDE  39  19  53.00  LONGITUDE  77  11  35.00  LEVEL  VALUE  1.217 

4  78  7  6  78  7  6  WWV  4  1  28.000  4  1  28.000 

4  78  7  6  78  7  6  WWV  6  15  49.000  6  15  49.000 

CHRONOMETER  OFFSET  0  0  0.000 

CHRONOMETER  RATE  0.00000  SECONDS/MINUTE 


1 

676 

NE 

12 

9 

54.6 

14.8 

14.6 

4 

8 

41.355 

15.90C 

.2 

29.4 

1 

676 

NW 

347 

50 

39.0 

14.0 

15.2 

4 

10 

41.126 

29.471 

-1.2 

29.2 

0 

684 

NW 

357 

10 

36.8 

14.3 

15.0 

4 

27 

38.286 

-  0.00 

-  .7 

29.3 

2 

684 

NE 

2 

49 

23.4 

14.8 

14.6 

4 

29 

28.846 

-  0.00 

.2 

29.4 

3 

695 

NE 

33 

23 

7.8 

14.9 

14.4 

4 

33 

55.076 

15.60C 

.5 

29.3 

3 

695 

NW 

326 

36 

50.0 

14.0 

15.3 

4 

36 

8.696 

29.471 

-1.3 

29.3 

4 

1483 

NW 

353 

7 

41.1 

14.8 

14.6 

4 

45 

50.728 

-  0.00 

.2 

29.4 

4 

1483 

NE 

6 

52 

19.5 

14.9 

14.5 

4 

47 

36.161 

-  0.00 

.4 

29.4 

5 

1488 

SE 

347 

18 

47.2 

15.2 

14.2 

4 

57 

56.713 

-  0.00 

1.0 

29.4 

5 

1488 

SW 

12 

41 

10.6 

14.0 

15.5 

4 

59 

28.983 

-  0.00 

-  1.5 

29.5 

6 

705 

SW 

5 

59 

21.2 

15.7 

13.8 

5 

1 

56.975 

15.30C 

1.9 

29.5 

6 

705 

SE 

354 

0 

37.6 

15.0 

14.4 

5 

3 

40.786 

29.471 

.6 

29.4 

7 

709 

SE 

324 

51 

26.6 

14.5 

15.0 

5 

7 

56.375 

-  0.00 

-  .5 

29.5 

7 

709 

SW 

35 

8 

32.0 

14.7 

14.9 

5 

9 

24.263 

-  0.00 

-  .2 

29.6 

8 

719 

SW 

3 

15 

43.1 

14.7 

15.0 

5 

19 

7.406 

-  0.00 

-  .3 

29.7 

8 

719 

SE 

356 

44 

18.0 

14.0 

15.8 

5 

20 

48.604 

-  0.00 

-  1.8 

29.8 

9 

723 

NE 

28 

17 

9.4 

14.9 

14.7 

5 

25 

6.539 

-  0.00 

.2 

29.6 

9 

723 

NW 

331 

42 

48.0 

14.2 

15.2 

5 

26 

52.965 

-  0.00 

-  1.0 

29.4 

10 

729 

NW 

326 

1 

25.6 

14.9 

14.8 

5 

28 

29.914 

0.00 

.1 

29.7 

10 

729 

NE 

33 

58 

33.7 

14.7 

15.0 

5 

30 

19.911 

-  0.00 

-  .3 

29.7 

11 

1506 

NE 

5 

33 

24.7 

15.2 

14.4 

5 

37 

43,384 

1 5.00C 

.8 

29.6 

11 

1506 

NW 

354 

26 

33.5 

14.7 

15.0 

5 

39 

38.108 

29.461 

-  .3 

29.7 

12 

1510 

SW 

4 

55 

15.5 

14.8 

15.0 

5 

43 

30.488 

-  0.00 

-  .2 

29.8 

12 

1510 

SE 

355 

4 

43.4 

15.5 

14.2 

5 

45 

10.503 

-  0.00 

1.3 

29.7 

13 

738 

NE 

10 

50 

17.0 

15.0 

14.6 

5 

48 

26.230 

-  0.00 

.4 

29.6 

13 

738 

NW 

349 

9 

41.0 

13.6 

16.3 

5 

50 

1.141 

-  0.00 

-2.7 

29.9 

14 

741 

SW 

28 

45 

37.8 

14.6 

15.2 

5 

57 

52.125 

-  0.00 

-  .6 

29.8 

14 

741 

SE 

331 

14 

21.2 

16.2 

13.3 

5 

59 

15.669 

-  0.00 

2.9 

29.5 

15 

749 

SE 

327 

1 

59.3 

14.9 

14.8 

6 

6 

52.920 

-  0.00 

.1 

29.7 

15 

749 

SW 

32 

57 

58.2 

14.0 

15.7 

6 

8 

14.614 

-  0.00 

-1.7 

29.7 

16 

1523 

SW 

11 

37 

59.4 

14.1 

15.6 

6 

12 

43.687 

15.00C 

-  1.5 

29.7 

16 

1523 

SE 

348 

22 

2.6 

14.3 

15.6 

6 

14 

16.649 

29.461 

-  1.3 

29.9 

SOURCE:  NOAA 
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Star 


Reduction  of  Zenith  Distance  to  the  Meridian 

Hour  Refr.  d  H  A7  n  — i  ■» *.• 


12°  09' 
347  SO 


5  7- ".88 
36  .12 


684 

•  357 
4/ 

10 

49 

i  36  .8 
23  .4 

-55  .35 
55  .51 

-  2  .80 
2  .80 

-  .43  ".84 

•12  -  .83 

357 

2 

10 

49 

I  34  .41 
’  25  .43 

695 

33 

326 

23 

36 

07  .8 
50  .0 

-68  .36 
65  .63 

36  .98 
-36  .98 

.30 
-  .79 

.72 

.67 

33 

326 

23 

36 

44  .36 
12  .90 

1483 

353 

6 

7 

52 

41  .10-52  .69 
19  .50  53  .03 

6  .76 

6  .76 

.12 

.24 

.77 

.76 

353 

6. 

7 

52 

35  .23 
25  .74 

1488 

347 

12 

18 

41 

47  .2 
10  .6 

-45  .54 
46  .98 

-12  .62 
12  .62 

.61 

-.91 

-.48 

.45 

347 

12 

18 

41 

34  .81 
22  .76 

705 

5 

354 

59 

00 

21  .2 
37  .6 

-49  .42 
54  .68 

5  .89 
-  5  .89 

1.16 
.37  - 

.61 

.75 

5 

354 

59 

00 

28  .86 
31  .33 

709 

324 

35 

51 

8 

26  .6 
32  ,0 

-40  .62 
47  .51 

-39  .48 
39  .48 

-  .30  - 

-  .12 

.09 

.07 

324 

35 

50 

9 

46  .73 
11  .43 

719 

3 

356 

15 

44 

43  .1 

18  .0 

-51  .17 
50  .31  ■ 

3  .20 
-  3  .20 

-  .18 
-1.10  - 

.68 

.68 

3 

356 

15 

44 

46  .80 
13  .02 

723 

28 
33 1 

17 

42 

9  A 
48  .0 

-54  .33 
52  .38 

30.18 

30.18 

.12  - 
-  .61 

.57 

.63 

28 

331 

17 

42 

9  .40 
48  .00 

729 

326 

33 

1 

58 

25  .6 
33  .7 

-56  .90 
53  .40 

-37  .80 
37.80  - 

.06 

-  .18  - 

.49 

.43 

326 

33 

00 

58 

48  .35 
10  .89 

1506 

5 

354 

33 

25 

24  ,7 

33  .5  ■ 

57  .63 
-57  ,41  - 

5  .46 

-  5  .46  - 

.49  - 
-  .18 

.91 

.90 

5 

354 

33 

26 

29  .74 
28  .75 

1510 

4 

355 

55 

4 

15  .50- 
43  .30 

-50  .36 
49  .93  - 

4.83  - 
-  4.83 

-  .12 
.79  - 

.64 

.63 

4 

355 

55 

4 

20  .85 

38  .75 

738 

10 

349 

50 

9 

17  .0  - 
41  .0 

-47  ,92  10,74 

47  .25  -10.74  - 

.24  - 
-1.64 

.62 

.60 

10 

349 

50 

9 

27  .56 

29  ,22 

741 

28 

331 

45 

14 

37  .8  - 
21  .2 

-41  .89 

41  .88  - 

30.78  - 
-30.78 

-  .37 

1.76  - 

.17 

.17 

28 

331 

46 

13 

08  .38 

52  .01 

749 

327 

32 

1 

57 

59  .3  - 
58  .2 

-40  .75  - 
41  .17 

•36  .37 

36  .37 

.06  - 
•1.03 

.09 

.09 

327 

32 

1 

57 

22  .89 

33  .63 

1523 

1 1 

348 

37 

22 

59  .4  - 
2  .6 

-46  .77 

46  .45  - 

1 1  .35  - 
It. 35  - 

.91 

.79  - 

.38 

.38 

11  . 
348  ; 

38 

21 

10  .22 

50  .08 

R  =  A0(n0  -  1)  tanZ  +  A,(n0  -  1)  tan3Z 


The  meanings  of  columns  in  table  1  are 


Column  1  = 

Column  2  = 
Column  3  = 
Column  4  = 
Column  5  = 
Column  6  = 
Column  7  = 
Column  8  = 


From  equation  (24) 

e  +  p  = 


Sequence  of  observations 

Star’s  number 

Observed  zenith  distances 

Level  readings,  in  the  order  left-right 

Universal  time  (WWV) 

Temperature  and  barometer  readings 
Difference  between  level  readings 
Sum  of  level  readings 

245°  32  '  33  ”.83 
245°  32  '  15  ".36 


-  -4~  ISZd  -  2Zr! 


e  +  p 


18  '47 
30 


0  ”.66 


Adding  the  value  (e  +  p)  to  the  Zd  and  subtracting  it  from  Zr,  we  get  the  values 
(Zd  +  c  +p)  and  Zf  -  (c  +  p)  and  are  shown  in  column  4  of  table  3.  By  subtracting 
these  values  (column  4)  from  each  mean  (column  5),  the  V  values  are  obtained. 
The  probable  error  of  an  observation  is  obtained  by  squaring  and  adding  the  V  values 
and  applying  equation  (17). 

/ IV  V| 

e  =  p  y/ - 

n 


The  sum  of  the  V  values  squared  gives 

|V  V]  =  26.7781 


and 


TABLE  3.  Determination  of  Observation  Error 


Star  Zenith  Dist  2(e+p)+V,  +V2  Zd+e+p  Means  V 

Zd/Zr  Zr~(e+p) 


684 

2° 

49  '25  ‘ 

'.49 

-0.10 

2° 

.49' 

24 

".83 

25.54 

+  .71 

2 

49 

25 

.59 

26 

.25 

-  .71 

695 

33 

23 

47 

.10 

33 

23 

46 

.44 

-  .71 

33 

23 

44 

.36 

+2.74 

45 

.02 

45.73 

+  .71 

1483 

6 

52 

25 

.74 

+0.97 

6 

52 

25 

.08 

25.26 

+  .18 

6 

52 

24 

.77 

25 

.43 

-  .17 

1488 

12 

41 

22 

.76 

-4.43 

12 

41 

22 

.10 

23.98 

+  1.88 

12 

41 

25 

.19 

25 

.85 

-1.87 

705 

5 

59 

28 

.86 

+0.09 

5 

59 

28 

.20 

28.76 

+  .56 

5 

59 

28 

.67 

29 

.33 

-  .57 

709 

35 

09 

1 1 

.43 

-1.84 

35 

09 

10 

.77 

12.35 

+  1.58 

35 

09 

13 

.27 

13 

.93 

-1.58 

719 

3 

15 

46 

.80 

-0.18 

3 

IS 

46 

.14 

46.89 

+  .75 

3 

15 

46 

.90 

47 

.64 

-  .75 

723 

28 

17 

39 

.13 

-3.03 

28 

17 

38 

.47 

40.64 

-2.17 

28 

17 

42 

.16 

42 

.82 

+2.17 

729 

33 

58 

10 

.89 

-0  76 

33 

58 

10 

.23 

11.27 

1.04 

33 

58 

I  1 

.65 

12 

.31 

-1.04 

1 506 

5 

33 

29 

.74 

-1.50 

5 

33 

29 

.08 

30.49 

1.41 

5 

33 

31 

.24 

31 

.90 

-1.41 

1510 

4 

55 

20 

.85 

-0.42 

4 

55 

20 

.19 

21.06 

0.87 

4 

55 

21. 

27 

21 

.93 

-0.87 

738 

10 

50 

27 

.36 

-3.42 

10 

50 

26 

70 

29.07 

2.03 

10 

50 

30 

.78 

31 

,44 

-2.03 

741 

28 

46 

08 

.38 

+0.39 

28 

46 

7 

.72 

8.18 

0.46 

28 

46 

07 

.99 

8 

.65 

-0.46 

749 

32 

58 

33 

.63 

-3.48 

32 

58 

32 

.97 

35.37 

2.40 

32 

58 

37 

.1  1 

37 

.77 

-2.40 

1523 

1  1 

38 

10 

->  i 

+0.30 

1 1 

38 

9 

.60 

10.09 

0.49 

1  1 

38 

09 

7)2 

10 

.58 

-0.49 

The  probable  error  of  an  observation  is 

e  =  0.6745  ^26.7781/15  =  ±  1  ''.20 


Further  investigations  using  random  stations  show  large  discrepancies  owing 
to  changes  in  the  TALCOTT  level  position  from  one  star  to  the  next.  However, 
these  discrepancies  do  not  affect  the  mean  value  of  the  individual  star  zenith  distances 
that  result  from  the  direct  and  reverse  instrument  position  because  of  the  cancellation 
in  the  mean  value.  As  a  result,  the  probable  error  of  an  observation  cannot  be  used 
by  this  method.  Thus,  another  procedure  must  be  used.  One  such  procedure  is  dis¬ 
cussed  next. 


Probable  Error  of  One  Observation  as  Function  of  Zenith  Distances. 

The  probable  error  can  be  computed  by  computing  the  zenith  distances  of  any  north¬ 
ern  star  with  all  southern  stars,  and  any  southern  star  with  all  of  the  northern  stars. 
The  zenith  distances  are  reduced  to  the  meridian  plane  as  shown  previously.  The 
condition  becomes 


n 


=  6 


1  n 


-  K  - 


"it 


for  the  northern  star,  and 


Zl*  -  Sli  “  5in  ~  Zin 


for  the  southern  star. 


Note  that  only  two  stars  can  be  chosen.  The  V  values  are  formed  with 
respect  to  the  means  of  each  star  (see  table  4). 
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$ 
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n  -  m  =  14  IV  VI  =  5.0137 

l. 

Error  of  one  observation 


INFLUENCE  OF  REFRACTION  IN  MEAN  LATITUDE 
VALUE  BY  THE  STERNECK  METHOD 


In  the  measurements  of  astronomic  coordinates,  a  degree  of  uncertainties  is 
due  to  a  lack  of  knowledge  of  the  atmospheric  refraction  during  the  acquisition  of 
stellar  positions.  Atmospheric  refraction  comprises  both  astronomic  and  terrestrial 
refraction;  however,  only  the  effect  of  astronomic  refraction  is  considered.  Astro¬ 
nomic  refraction  is  defined  as  the  bending  of  light  from  objects  at  distances  above  the 
earth’s  mean  radius.  From  a  star,  the  light  moves  as  it  passes  through  strata  of  different 
densities.  Turbulence  in  the  earth’s  atmosphere  causes  changes  in  the  brightness  of  a 
star.  Small  angular  displacement  can  be  expected  during  optical  observations. 


Since  the  determination  of  the  actual  air  density  profile  is  difficult,  it  is  gener¬ 
ally  approximated  by  using  a  quasi-standard  model.  The  solution  of  the  refraction 
integral  may  then  be  represented  in  the  form 

R  =  A  tanZ  +  B  tan3Z  +  CtansZ  +  Em  (25) 

where  the  coefficients  A,  B,  and  C  depend  on  the  quasi-standard  model.  There 
remains,  however,  an  error  owing  to  the  replacement  of  the  actual  air  density  profile 
by  a  model  profile,  Using  the  former  requires  practically  the  numerical  evaluation  of 
the  refraction  integral. 


In  geodetic  astronomy,  the  first  two  terms  are  considered  that  give  satisfactory 
results  to  zenith  distances  up  to  Z  =  75°. 

The  stars  used  for  latitude  determination  are  chosen  in  a  way  with  no  zenith 
angle  greater  than  30°.  For  the  mean  latitude  value  that  results  from  a  night  of 
observation,  the  latitude  is 

*  =  —  (26n  +  25s)  +  (2Zn  -  2Z.)  +  dip  (26) 

n  n  n  * 

where 

dp  =  — A.  (£  tan  Zn  -  2  tan  Zs)  +  -A-  (2  tan3  Zn  -  2  tan3  Zs)  (27) 
n  n 

and  n  represents  the  number  of  observed  stars.  A  program  for  latitude  comprises 
16  stars,  which  are  evenly  divided  to  north  and  south  of  the  zenith.  The  coefficient 
B  is  small,  less  than  0.07  ".  Since  no  zenith  distance  is  greater  than  30°,  the  second 
term  of  equation  (27)  vanishes  for  practical  purposes. 
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ERROR  LEVEL  CORRECTION  •  A  TALCOTT  level  is  used  by  NOAA 
for  zenith  distances  observations.  The  graduation  readings  increase  outwardly  from 
each  side  from  the  center  zero  mark  on  the  .vial-level  scale. 

Figure  3  shows  the  position  of  the  instrument  with  an  inclination  i. 


FIGURE  3.  The  Inclination  i. 

The  bubble  reads  higher  to  the  right  and  lower  to  the  left.  Let  R  and  L  equal  the 
right  and  left  readings  and  the  mean  reading  is  denoted  by  L.  Let  p  be  the  level 
sensibility. 


The  inclination  i  is 

i  =  p(R  -  L0 )  =  p(L  -  L„)  (28) 

The  direct  zenith  distance  for  0°  to  90°  is  greater  by  i.  and  the  reverse  Zr 
for  0°  to  270°  is  lesser  by  i.  For  Zj,  the  correction  for  inclination  is  negative. 
Then.  R  is  greater  than  L0.  therefore 

i  =  p(L()  -  R)  =  -  (L  -  R) 
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where  the  direct  zenith  distance  is 

i  =  4  (L-R) 

tm. 

and  the  reverse  zenith  distance  is 


i  =-y  (L-R) 


When  the  observer  maintains  a  small  difference  (L  -  R),  an  error  in  the  value 
of  p  has  little  influence.  For  an  error  5p  =  0".l  and  (L-R)  =  3,  the  error  in 
the  inclination  is  dj  =  0  ” .  1 5. 


CONDITIONAL  EQUATION  •  When  summing  equations  (18)  and  (19), 
the  mean  value  is  free  of  the  index  error  and  personal  equation,  as  shown  in  equation 
(20). 

Z  =  —  (Zd  +Zr)  +  i-(V,  -V2) 


Let  index  n  refer  to  the  north  star  and  index  s  refer  to  the  southern  star. 
The  sum  of  all  of  the  star’s  zenith  distances  have  to  satisfy  the  equation 

IZ„  +  £ZS  +  SV,  -  SV2  =  S5n-£5s  (29) 

Station  values  used  to  determine  the  observational  error  are  applied  to  equation 
(29)  as  follows: 


Z6„ 

L5S 


-  -  448° 

35 

'  " 

.53 

*4 . 

-  -133° 

56  ’ 

17. 

r  i 
O 

-  •  1  79 

15 

02 

.09 

*4 . 

-  -135 

24 

0b. 

79 

=  269° 

20 

'20" 

.44 

IZn  +  2ZS 

=  269° 

20' 

23. 

”81 

Then 

SV,  -  IV2  =  -  3  ".37 
The  sum  of  V(  and  V2  is  inseparable. 
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LEAST  SQUARES  SOLUTION  •  The  discrepancy  shown  in  equation 
(26)  can  be  minimized  if  a  solution  by  a  least-square  method  is  followed.  Two  sets 
of  equations  of  condition  are  established  that  are  based  on  previously  derived  general 
equations  having  all  stars  observed,  or  reduced,  to  a  constant  hour  angle.  The  equation 
is 


sin  ip  sin(5n  -  6S)  +  cos  6n  cos  Zs  =  cos  5S  cos  Zn 


(30) 


The  index  n  and  s  refer  to  the  north  and  south  star,  respectively.  Since  the 
zenith  angles  are  small,  the  formula  fulfills  the  accuracy  requirement.  Assuming  small 
errors  in  Zs  and  Zn  and  differentiating  relative  to  0  and  Zs,  Zn  yields 


COS  ip  dip 


cos  6n  sin  Zs 
sin  (6n  -  5S) 


cos  sin  Zn 

S  It  4  rj 

sin  (5n-  5S) 


(31) 


Since  the  sin  Z’s  are  always  small,  the  errors  in  the  zenith  distances  decreases  with 
th.e  factor  sin  Z.dZ. 


The  least-square  method  continues  as  follows: 

A  star  with  large  declination  is  selected  from  the  north  set.  The  zenith  distance  for 
the  star  is  unknown.  The  star  is  used  with  each  of  the  stars  from  the  south  to  obtain 
a  number  of  free  equations  of  the  form 

\  -  X  +  a,  Y,  =  «,  (32) 

k  where  the  unknowns  X  and  Y  ]  are 

X  =  sin  ip 

Y,  =  cos  Zn  (33) 


and  the  coefficients  ajt  Kj,  are 


Similarly,  by  holding  a  fixed  southern  star  with  each  one  from  the  north  set. 
one  can  obtain  m  observation  equations  such  as 

X  +  b.  Y2  =  £.  (36) 

where  the  unknown  Y2  is 

Y2  =  cos  Zs  (37) 

and  the  coefficients  b’s  and  id’s  are  obtained  from 


cos  6in 
sin(6jn  -6S) 


cos  6S  cos  Zin 
sin  (6in  -  6S) 


(38) 


j  -  1.2 . ,m 

Solving  equations  (32)  and  (36)  simultaneously  by  the  method  of  least 
squares  and  eliminating  the  unknown  X,  one  obtains  two  normal  equations  as 

AY,  +  BY2  =  L, 

BY,  +  CY2  =  L2  (39) 

To  these  normal  equations  is  added 

Z,  +  Z2  =  6n  _  6S  (40) 

L.et  Z,  and  Z2  be  the  observed  values  and  dZ,.  dZ2  be  the  corrections 
needed  to  satisfy  equation  (40),  and  let 

WI2  =  5n  “  6s  ~  (ZI  +Z2>  (41) 

then 

dZ,  -!  dZ2  =  W,2  (42) 


To  solve  equation  (40)  with  equation  (42),  proceed  as  follows. 


By  replacing  the  unknowns  Y,,  Y2  by 

Yj  =  cosZj  -dZj  sinZ, 


r  I  =  eos£[  —  uz-j  sin  £.| 

f2  -  cos  Z2  -  dZ2  sin  Z2 


equation  (40)  becomes 

A  sin  Z,  dZj  +  B  sin  Z2  dZ2 
B  sin  Zj  dZj  +  C  sin  Z2  dZ2 


=  W, 


where 


=  (A  cos  Zj  +  B  cos  Z2  -  Lj  )/sin  1 
=  (B  cos  Zt  +  C  cos  Z2  -  L2)/sin  1 


The  normal  equation  (44)  is  to  be  solved  by  using  equation  (43). 

It  is  solved  by  introducing  the  unknown  K  as  follows: 

A  sin  Zt  dZj  +  B  sin  Z2  dZ2  +  K  =  W, 

B  sin  Z,  dZ[  +  C  sin  Z2  dZ2  +  K  =  W2 
dZ,  +  dZ2  =  W12 


Solving  for  dZ(  and  dZ2  and  using 
Yj  =  cos (Zj  +  dZj ) 


Y,  =  cos (Z2  +  dZ,) 


The  unknown  X  is  computed  from  the  first  normal  equation,  and  the  latitude  is 
computed  from 

\p  =  arc  sin  X 

The  test  conducted  refers  to  the  Station  Aero  USAETL  1978.  Star  number 
676  is  chosen  for  the  north  star  and  star  1488,  for  the  southern  set.  The  data  used 
is  shown  in  table  5. 
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TABLE  5.  Final  Meridian  Zenith  Distances 


STAR 


1 

676 

17 

56 

5.11 

17 

56 

8.43 

51 

29 

42.20 

12 

9  49.30 

\ 

W 

1 

684 

18 

15 

.54 

18 

15 

.46 

42 

9 

18.06 

2 

49  25.49 

© 

3 

695 

18 

21 

29.92 

18 

21 

31.29 

72 

43 

38.94 

33 

23  45.61 

1 

4 

1483 

18 

33 

13.40 

18 

33 

13.23 

46 

12 

18.18 

6 

52  25.24 

■lit 

5 

1488 

18 

45 

14.77 

18 

45 

14.05 

26 

38 

29.45 

12 

41  24.00 

m 

6 

705 

18 

49 

21.48 

18 

49 

18.85 

33 

20 

24.48 

5 

59  28.75 

■M. 

7 

709 

18 

55 

13.88 

18 

55 

10.43 

4 

10 

40.14 

35 

9  12.27 

■i 

*K»s 

8 

719 

19 

6 

33.42 

19 

6 

33.85 

36 

4 

7.13 

3 

15  46.87 

f 

9 

723 

19 

12 

36.16 

19 

12 

37.13 

67 

37 

33.96 

28 

17  40.64 

$ 

10 

729 

19 

16 

1.88 

19 

16 

3.63 

73 

19 

5.46 

33 

59  11.20 

f 

11 

1506 

19 

25 

19.24 

19 

25 

19.35 

44 

53 

23.78 

5 

33  30.48 

1 

12 

1510 

19 

30 

59.92 

19 

31 

.13 

34 

24 

32.48 

4 

55  21.06 

13 

738 

19 

35 

53.91 

19 

35 

54.24 

50 

10 

21.95 

10 

50  29.06 

14 

741 

19 

45 

15.66 

19 

45 

15.66 

10 

33 

46.19 

28 

46  8.15 

a 

,4 

15 

749 

19 

54 

17.00 

19 

54 

16.79 

6 

21 

19.00 

32 

58  35.33 

16 

1523 

20 

0 

14.38 

20 

0 

14.54 

27 

41 

44.22 

11 

38  10.36 

n 

The  following  equations  of  observation  are 
-X  +2.1 2664562  Y, 


+  2.1 2664562  Yt 
+  2.68118197 

=  1.44509638 
=  1.98718370 

I 

+  1.35671000 

=  .69244814 

+  3.03876241 

=  2.33673594 

Vcr 

+  2.80804611 

=  2.11119945 

! 

+  1 .50046665 

=  .83297305 

i 

+  1.40210030 

=  .73681725 

+2.19417462Y! 

=  1.51  110765 

1 

+  1.481 27975  Y2 

=  2.07890158 

% 

i 

+  2.77169144 

=  3.33779276 

+  .41216626 

=  1.03590425 

i 

-A 

+  2.06682246 

=  2.65014159 

1 

+  .53038532 

=  1.20001419 

l 

+  .39458299 

=  1.01875357 

§ 

i 

+  2.26246721 

=  2.84100790 

+  1.60419822Y2 

=  2.19881757 

t 
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Normal  Equations 

16X  -17.10808768Y,  +  11.57359364Y2  =  4.70777185 
39.70290799Y,  =  27.96836505 

22.50285187Y2  =  29.28858325 

Eliminating  the  unknown  X  yields 

2 1.409991 49 Yj  +  12.375 12842Y2  =  33.00217590 
1 2.375 1 2842  Y  2  +  14. 1 310975 1 Y2  =  25.88321835 


From  these  equations  and  from  the  data  in  table  5,  the  result  is 

W12  =  0  ".55 

W(  =  1  .36 

W2  =  -0  .14 

From  equation  (45),  the  result  is 

4.511  dZ,  +  2.7 19  dZ2  +  K  =  -  1.36 

2.608  dZ,  +  3.104  dZ2  +  K  =  -0.14 

dZ,  +  dZ2  =  -0.55 

Then,  solving  for  dZt  and  dZ2. 

dZ,  =  -  0  ".687  Z,  =  Z,  +  dZ,  =  12°  09'  48  ".61 

dZ2  =  +  0  .137  Z2  =  Z2  +  dZ2  =  12  41  24  .14 

With  these  values  of  Z,  and  Z2,  the  values  of  Y|  and  Y2  are  obtained.  From  the 
first  normal  equation  is  obtained 

X  =  sin  y  =  0.633800660 

y  =  39°  19'  53 ''.50 


i 

% 


l 


The  mean  latitude  from  Sterneck’s  method  for  16  stars  from  table  6  is 


=  39°  19'  53  ".40 

More  tests  must  be  done  using  different  star  pairs.  One  example,  which  follows,  is 
observing  several  zenith  distances  in  a  direct  and  reverse  theodolite  position. 


A  star  whose  declination  is  defined  by  the  equation 

t  —  /^r  sin  <45  +  W 

tan  6  >  */  2  - - 

cos  <f> 

for  a  star  north  of  the  zenith 

♦  s  —  nr  sin  <45  ~  0) 
tan  6  <  v  2  - - 

cos  <t> 

for  a  star  south  of  the  zenith,  has  a  very  small  variation  in  zenith  distance  with  res¬ 
pect  to  the  minimum  zenith  distance  when  it  crosses  the  upper  meridian.  Hence, 
if  we  observe  it  continuously  and  record  the  zenith  distance  with  the  corresponding 
times  for  approximately  10  to  15  minutes  before  and  after  it  crosses  the  meridian 
and  continue  to  observe  it  for  about  the  same  time  after  it  has  crossed  the  meridian, 
the  star’s  image  will  always  be  in  the  field  of  view  of  the  telescope.  Thus,  many  ob¬ 
servations  can  be  made  in  a  direct  and  reverse  theodolite  position  by  using  the  tangent 
screw  of  the  vertical  circle  or  by  measuring  the  small  angles  with  the  eyepiece  micro¬ 
meter.  These  zenith  distances  can  be  reduced  to  the  meridian  by  using  equation  (2.1) 
to  (2.6).  The  average  value  is  then  free  from  index  error  and  personal  equation  error. 
The  refraction  corrections  are  practically  the  same  for  all  observations  since  the 
changes  in  zenith  distances  reach,  at  most,  3  to  5  minutes. 


MfSr 

l 


I 
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TABLE  6.  Results  of  Latitude  Determined  by  the  Sterneck  Method 

1  STATION  AERO  USAETL  1978  MD  INST  NO  T-4  86994 
5  «5  «  «•  ACC  NO  NE°CCT  4'Tw 


4 

5 


78  7  6  78 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
l(> 


7  6 

wwv 

4 

1 

28.000 

4 

1  28.000 

*  T  /I 

7  6 

wwv 

6 

15 

49.000 

6 

15  49.000 

676 

N 

39" 

19' 

52  ".90 

.50 

684 

N 

39 

19 

52  .57 

.83 

69  S 

N 

39 

19 

53  .34 

.06 

1483 

N 

39 

19 

52  .94 

.46 

1488 

S 

39 

19 

53  .45 

-  .05 

705 

S 

39 

19 

53  .23 

.17 

709 

s 

39 

19 

52  .41 

.99 

719 

s 

39 

19 

54  .00 

-  .60 

723 

N 

39 

19 

53  .32 

.08 

729 

N 

39 

19 

54  .27 

-  .87 

1 506 

N 

39 

19 

53  .30 

.10 

1510 

S 

39 

19 

53  .54 

-  .14 

738 

N 

39 

19 

52  .89 

.51 

741 

S 

39 

19 

54  .34 

-  .94 

749 

S 

39 

19 

53  .34 

.06 

1523 

s 

39 

19 

54  .58 

-1.18 

MEAN  LATITUDE  IS  39  19  53.400 
SID  ERR  SlNCiLE  OBS  IS  .6240  SECS 
STD  ERR  OE  Till;  MEAN  IS  .1500  SECS 
ECC  LAT  IS  0.00  SECS 

KCC  EON  IS  0.00  SECS 


USING  16  ACCEPTABLE  OBSERVATIONS  FROM  THE  16  OBSERVE 

STATI0N  ANU  A,’I>LV,NC  THE 


THI*  ADJUST lil)  LATITUDE  IS  39°  |9*  53  "A0Q 


t 


{ 
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EVALUATING  LATITUDE  FROM  A  SET  OF  NORTH  STARS  •  A 

test  was  conducted  for  Aero  Station  using  the  author's  formula 


sin  $ 


cos  Z]  cos  6 2  -  cos  Z2  cos  6] 
sin  (62  -  5( ) 


(46) 


Star  number  1  is  held  fixed  with  respect  to  the  remainder  of  the  same  set. 
For  star  number  2,  the  large  declination  is  used.  For  the  north  set,  star  number  676 
and  tor  the  south  set,  star  number  719  arc  selected.  Results  are  shown  in  table  7. 


TABLE  7.  Evaluating  Latitude  By  Set 
North  Set  South  Set 


mean 


39° 

19' 

53 

"43 

39° 

19  ' 

53 

".23 

53 

.36 

52 

.24 

53 

.37 

52 

.18 

53 

.35 

52 

.58 

53 

.33 

54 

.40 

39 

19 

53 

.41 

53 

.23 

39 

19 

54 

.84 

39 

19 

53 

.38 

39 

19 

53 

.24 

I 


t 
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The  results  of  these  two  sets  show  that  better  quality  of  observation  was 
achieved  from  the  north  stars.  The  differences  may  be  caused  by  sky  illumination 
and  refraction  anomalies.  More  observations  are  recommended  for  testing  the  different 
methods. 


STATION  LATITUDE  ERROR  •  The  generalized  form  of  Gauss’s  pro¬ 
pagation  error  law  is 


in 


*  *  (-sk  m')2  +  (irms) 


2  „  3f  3f 

+  2  -  -  m,  * 

3Z  36  10 


(47) 


where  m26.  the  covariance  of  Z  and  6,  is  given  by 

'»z3  “  P  •  n> z  '  m6 
Because  m2  and  mg  are  independent. 


,nz6  =  0 


In  the  Sterneck  latitude  method, 


3f  _  3f 
3Z  "  35 


therefore 


n»J  "  ml  +  ml 


(48) 


The  standard  error  of  the  FK4  system  of  the  FK4  for  declination  is  1,925; 
so  the  standard  error  in  declination  is  calculated  from  formula 


M  =  (ai\.  +  (^rsr1)2  < 


For  Aero  Station  T  =  1,958,  so  we  have 


("4, =  (™4 


+  0.281  m:. 

T  I* 

1  n 


For  the  mean  value  of  the  latitude  derived  from  1 6  stars, 

.016 


(■-). 


=  0  .06 


(“4 


=  0.04 


The  standard  error  of  an  observation  is  mz  =  0.60,  with  n  =  16.  the  result 


is  mz  =  0.15. 


31 


The  standard  error  of  the  mean  latitude  is 


±  0  ".16 


For  two  nights  at  each  station,  the  average  value  for  the  entire  test  course  was 
between  0.18  to  0.20.  The  most  frequent  value,  0.20,  is  used  as  the  mean  error 
of  latitude  for  any  given  station. 


IMPROVING  ACCURACY  IN  LATITUDE  DETERMINATION  •  Improving 
the  accuracy  in  latitude  can  be  achieved  by  reducing  all  stars  to  a  common  hour  angle, 
i.e.  t  =  60s.  The  east  and  west  reduced  zenith  angles,  having  the  mean  value  free 
of  the  index  error,  are  the  same.  With  respect  to  the  same  hour  angle,  one  star  from 
the  north  and  two  from  the  south  are  selected  with  indices  1,  2,  and  3.  Let 

6,  >  62  >  > 

An  equation  of  condition  that  is  free  of  the  latitude  and  is  related  to  the  zenith 
angles  and  declinations  is  derived  as  follows: 

cosZ,  sin(S2-83)  +  cosZ3sin(8,  -82)  *  eosZ2sin(8,  -63)  (49) 

with 


z. 

+  Z3 

=  8,-83 

+  8Z, 

+  6Z 

z. 

+  Z2 

=  8,  -  82 

+  6Z, 

+  8Z. 

Similarly,  selecting  two  stars  from  the  south  and  one  from  the  north  with  the 
condition 


6,  >  82  >  83  > 

and  with  the  equation  (49)  the  condition  is 

Z2  +  Z3  =  6Z  -83  +  6Z2  +  8Z3 

Zj  +  Z3  =  8j  -83  +  8Zj  +  8Z3 


(51) 
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Tlie  zenith  distances  in  equations  (30)  to  (33)  are  obtained  as  follows: 

With  the  original  hour  angles,  the  observed  zenith  distances  are  corrected  by  applying 
equations  (14)  and  (15) 

In  figure  4,  the  configuration  is  depicted  for  a  star  north  of  the  zenith  6  >  0, 
where  the  circle  SM  is  the  star  trajectory  and  S  is  the  star  position. 


NP 


Z 


FIGURE  4.  The  Configuration  for  a  Circumpolar  Star. 


Let  t  be  the  selected  common  hour  angle. 


ZM  =  Zn, ,  Meridian  Zenith  Distance 


zs  =  z, 

SP  =  MP  =  90°  -6V 


then 

from 


S'P  =  90°  -6, 


tan  6.  =  tan  6  sec  t 


(52) 


resulting  in 


>  v  >,  ■*  y  ^  y.>  y  -.vr-  -w  .*««*■'*'*-*'*-'*  |B  MMIHKMRRI^IKW |Ti!*fi(i(SNSj533^,’^^r^^r" 


cos  Zt 


cos  ZS 


sin  5 
sinSx 


(53) 


where 


ZS’ 


Z.„  +  6. 
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LONGITUDE  DETERMINATION  BY  TRACKING  STARS  CROSSING  THE 
MERIDIAN  PLANE.  •  Longitudes  were  determined  using  a  Wild  T  -  4  eo- 
dolite.  The  observer  tracks  the  star  halfway  toward  the  center  of  the  field  of  ew 
with  the  eye  piece  micrometer.  The  instrument  is  then  reversed  180°,  the  telescope 
reset  for  the  same  star,  and  the  observer  resumes  tracking  the  star.  This  procedure 
of  tracking  the  star  before  and  after  crossing  the  meridian  must  be  discussed. 


The  procedure  of  reversing  a  transit  instrument  by  lifting  it  from  its  bearings 
to  track  the  star  before  and  after  crossing  the  meridian  plane  is  different  from  the 
procedure  used  with  the  Wild  T  -  4  in  which  the  theodolite  is  rotated  along  a  vertical 
axis  180°  with  respect  to  the  first  position  through  the  horizontal  scale  readings. 
This  procedure  results  in  errors  which  influence  the  longitude  or  time  computations. 
Errors  requiring  analysis  are 

1 .  Error  of  a  transit  observation. 

f  2.  Instrumental  errors, 

i  3.  Rotation  of  the  pier,  if  any. 

\  4.  Equation  of  bright. 


I .  Error  of  a  Transit  Observation.  The  error  in  the  observed  time 
of  star  transit  is  independent  of  the  personal  equation  and  other  constant  errors.  The 
final  result  is  affected  by  the  error  when  the  observations  of  the  two  observers  are 
not  combined.  The  error  is  determined  by  comparing  the  values  of  the  thread  intervals 
or  micrometer  position  from  the  observations.  Let  r  equal  the  equational  interval 
between  two  closely  ruled  vertical  lines  on  glass,  and  let  Tj,  T2  equal  the  observed 
transit  times.  Then 


.  WJ 

••J 

-.4, 

•£jm : 

y 

•  -1 
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bach  observation  furnishes  a  value  of  r.  When  considering  a  great 
number  of  values,  the  error  of  a  single  determination  is  obtained  from 


/  tv  V)" 

V  (m  -  1 ) 


(55) 


in  which  the  values  of  V  are  the  residuals  computed  by  subtracting  the  known  value 
of  r  from  each  value  in  the  M  number  of  observations.  The  error  e  of  the  observed 
time  transit  over  each  thread  is 


nv 

if 

'I 

I 

f. 

t. 


and 


2e2  =  r2 


(V  VJ 


2(m  -  1) 


(56) 


If  n  indicates  the  number  of  thread  with  the  mean  for  a  star  having  declination 
6  =  0,  the  error  is 


£ 


[VV1 


2(m  -  l)n 


(57) 


For  a  star  with  nonzero  declination,  the  error  is 


€ 


(V  VI 

2(m  -  1  )n 


sec  6 


(58) 


2,  Instrumental  Errors  are  from  Defects  in  the  Instrument.  Also 
included  are  errors  made  in  setting  up  the  instrument  (eollimation,  inclination,  and 
azimuth  of  the  horizontal  axis)  and  imperfections  in  the  pivots.  A  good  observational 
procedure  can  reduce  the  instrumental  errors. 


3.  Rotation  of  the  Pier.  The  instrumental  errors  are  controlled 
by  placing  two  ground  marks  north  and  so-  ;  of  the  meridian  plane.  Any  changes 
will  affect  the  instrument  azimuth  orientc*io:i  and  influence  the  longitude  or  the 
time  determination. 


i 

t 


6 
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4.  Equation  of  Bright.  The  eye  perceives  the  light  of  a  dim  star 
later  than  that  of  a  bright  star.  To  overcome  the  difference  of  perception,  the  observer 
tracks  the  star  with  a  micrometer  that  has  two  parallel  threads  in  place  of  one.  The 
star  is  centered  in  the  middle  of  the  two  threads  and  then  tracked. 


Analysis  of  Data.  The  micrometer-contact  times  from  the  longitude 
observations  are  not  available  for  these  analyses.  Instead  NGS  (National  Geodetic 
Survey,  NOAA)  provided  a  computed  mean  time  for  a  star’s  meridian  transit.  The 
transit  time  error  is  estimated  in  this  report  by  methods  independent  of  instrument 
orientation  errors. 


The  star  is  tracked  before  and  after  crossing  the  meridian  by  revolving 
the  theodolite  through  the  vertical  axis.  The  observation  is  set  to  a  different  vertical 
plane  because  a  special  attachment  is  not  provided  to  the  theodolite. 


Before  the  star  meridian  transit,  the  time  refers  to  an  azimuth  a,  and  to 
an  azimuth  a2  after  the  star  crosses  the  meridian.  A  small  error  exists  in  the  meridian 
time  crossing. 


The  error  in  the  star  transit  time  is  given  by 


6T  = 


sin  (ip  -  5) 
cos  6 


(59) 


To  avoid  the  error,  stars  are  selected  between  declinations  that  are  close 
to  the  station  latitude.  The  effect  of  the  error  on  the  observed  transit  time  of  a  star 
set  is  examined.  The  longitude  is  determined  from  the  same  star  set. 


Inherent  Error  in  the  Transit  Time.  To  set  up  the  Wild  T  -  4  theodolite 
position  without  an  angular  deviation  when  reversing  to  observe  the  star  transit  time 
is  practically  impossible.  The  first  longitude  determination  of  Aero  Station  is  con¬ 
sidered.  The  deviation  of  azimuth  coefficient  results  from  a  fixed  north  star  related 
to  stars  from  the  south  as  follows: 

Pi  =  -  5X  +  a  •  Aj 

Pi  =  -UTi  +  0O  +  Adi  +  bB|  +  kt  -  (“j  -  X0 )  (60) 


36 


where 


SA  =  Correction  to  the  assumed  longitude  A0 
UTj  =  WWV  observation  time 
0O  =  Greenwich  Vernal  Equinox  hour  angle 
A  =  sin  —6)  seo  8 

B  =  cos  (i p  -  5)  see  6 

b  =  Level  error 

k  =  Correction  for  diurnal  aberration 
a  =  Azimuth  error  expressed  in  time 

■I-. 

•| 


Table  8  shows  the  data  used  in  this  investigation  and  table  9 
the  evaluation  of  mean  azimuth  orientation. 


Equations  to  be  used  are 
0i  =  -A5  +  a  Aj 


Pi  = 

U.Tj  +  0 

0  +  A02 

+  b  Bj 

+  kj  - 

(«i  +Xo) 

TABLE 

8.  Evaluation  of 

U.T 

“i 

8 

bB+K 

1 

23h 

'40m 

43s. 663 

21,'43m 

8s. 847 

9C 

'  46 ’.9 

+.053 

1 

43 

35  .064 

46 

1  .294 

49 

13  .1 

+.194 

3 

49 

39  .466 

21  52 

6  .212 

25 

49  .8 

-.007 

4 

23 

57 

35  .429 

22  0 

3  .650 

13 

1  .3 

+.235 

5 

0 

3 

32  .706 

6 

1  .874 

25 

14  .8 

+.133 

6 

7 

37  .684 

10 

8  .405 

58 

6  .2 

+.451 

7 

20 

12  .117 

22 

44  .682 

52 

7.7 

+.398 

shows 


Pi 


2s. 834 
3  .270 
2  .989 
2  .919 

2  .990 

3  .554 
3  .386 
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0O  =  3h07">  18s  .368 

X  =  5  8  49  .467 

A 0  =  (1  +  CQ)  U.T 
1+C0  =  1.0027379051 


The  coefficients  A  are 


Star 

A 

1 

.500 

2 

-  .263 

3 

.259 

4 

.455 

5 

.269 

6 

-.609 

7 

-.361 

TABLE  9.  Evaluation 

North 

South  \  2 

of  Azimuth 

6 

7 

1  -.571 

-.643 

-.640 

3  -.538 

-.651 

-.640 

4  -.489 

-.591 

-.572 

5  -.526 

-.642 

-.629 

Table  10  shows  that  star  numbers  2  and  4  have  errors  in  the  time 
crossing  because  of  misorientation.  A  graphic  representation  of  equation  (47)  shows 
a  deviation  of  the  stars  2  and  4  with  the  other  set  (figure  5). 


38 


where 


The  resulting  deviation  in  the  correction  X6  from  the  author’s  equation  is 


AT 


-  AX 


h  ’Ml  ~  ft/*  M) 

M,  -  M2 


(61) 


Mj  =  tan  $  -  tan  6, 
M2  =  tan  i p  -  tan  62 


(62) 


Solving  by  pairs  of  stars,  one  can  see  in  table  10  that  the  star  numbers 
2  and  4  have  an  erroneous  time  of  crossing,  which  is  caused  mainly  by  the  difference 
in  azimuth  orientation. 


TABLE  10.  Evaluation  Accuracy  of  AX 


Star 

1 

2 

3 

4 

5 

6 

7 

1 

3.120 

3.156 

3.171 

3.159 

3.155 

2 

3.120 

3.129 

3.132 

3 

3.15b 

3.081 

3. 1 58 

3.155 

4 

3.141 

3.081 

3.093 

3.191 

3. 1 79 

5 

3.171 

3.132 

3.093 

3.163 

3. 1 59 

6 

3.159 

3.159 

3.191 

7 

3.155 

3.155 

3.179 

I 

|\ 

1  I'hc  author  recommends  that  future  astro  teams  observe  pairs  of  stars  with 

|  instrument  position  clamp  east  (or  west)  and  pairs  of  stars  clamp  west  (or  east).2 

|  The  author’s  formulas  are  suggested  since  the  formulas  are  independent  of  instrumental 

a/imuth  orientation  and  reduce  the  number  of  unknowns  by  50  percent. 


|  _ 

I  ^National  Ocean  Survey,  formerly  the  (nasi  and  (leodctic  Surrey,  “Manual  oltieodetic  Astronomy,"  1947. 

t  Special  Publication  No.  2.17,  p.  54, 
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TABLE  11.  Results  From  The  Baidini  Method 
NS  S  N  N  S  N 


STAR 

v  Andr. 

a  Triang. 

a  Arietis 

4>  Persei 

5  Andr. 

7  Triang. 

r  Persei 

v  Antlromedae  N 

-0.729 

-0.720 

-0.736 

a  Trianguti 

S 

-0.736 

-0.738 

-0.743 

-0.736 

a  Arictis 

S 

-0.720 

-0.725 

-0.737 

-0.721 

4>  Persei 

N 

-0.738 

-0.725 

-0.743 

51  Antlromedae  N 

-0.744 

-0.737 

-0.745 

7  Trianguli 

S 

-0.736 

-0.743 

-0.745 

-0.742 

t  Persei 

N 

-0.736 

-0.721 

1 

-0.742 

MEAN 

-0.731 

-0.737 

-0.726 

-0.735 

-0.742 

-0.741 

-0.733 

AT  -  -  9s. 00  +  6T 
AT  =  -  9s. 735 


SOURCE:  National  Ocean  Survey,  formerly  the  Coast  and  Geodetic  Survey,  “Manual  of  Geodetic 

Astronomy,  1947,  Speciul  Publication  No.  237,  pp.  44  and  45. 

NOTE:  As  the  star  is  observed  before  and  after  crossing  the  meridian  plane,  the  effect  of 

coilimation  is  cancelled  out,  so  Buldini’s  method  solves  for  only  one  unknown:  AT. 
Each  star  can  be  combined  with  another  star,  forming  a  pair.  The  results  for  any  pair 
can  be  seen  at  intersecting  a  column  and  a  line  on  a  table.  The  values  are  the  decimal 
part  of  seconds  in  longitude.  The  final  value  is: 


AT  =  -A*  =  -9s  +  St  =  AT  =  -  9s. 735 


% 
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DEFLECTIONS  OF  THE  VERTICAL  •  The  £  and  n  components 
of  the  vertical  deflection  in  the  meridian  and  in  the  prime  vertical  are  related  to  the 
astronomic  latitude  and  longitude  <p>  X.  with  respect  to  the  geodetic  B  and  L, 
by  the  equations 


€  =  *  -  B 

n  =  (X  -  L)  cos  y 


(63) 


Table  1 2  gives  the  results  of  the  final  astronomic  positions  and  the  mean  errors 
in  its  positions.  Table  13  gives  the  errors  in  computing  £  and  n. 


TABLE  12.  Final  Astronomic  Positions 


STATION 

<i>  (N) 

°0 

A(W) 

°A 

AERO  USAETL  1978 

39° 

19' 

52 

".82 

±0 

”.17 

77° 

11  ' 

31 

".08 

±0 

".38 

BRINK  1969  AZ.  MK. 

39 

12 

27 

.98 

±o 

.19 

77 

14 

32 

.76 

±0 

.38 

CEDAR  HEIGHTS  1969 

39 

15 

14 

.73 

+0 

.20 

77 

13 

45 

.65 

±0 

.38 

CHEVY  USAETL  1978 

39 

04 

54 

.05 

to 

.20 

76 

56 

45 

.35 

to 

.38 

DALE  1943  RM  1 

39 

17 

21 

.22 

±0 

.20 

77 

15 

49 

.02 

to 

.38 

EIRE  USAETL  1978 

39 

17 

43 

.40 

±0 

.20 

77 

19 

00 

.77 

±0 

.38 

FREEWAY  USAETL  1978 

39 

1 1 

29 

.16 

±0 

.20 

77 

15 

04 

.02 

to 

.38 

HORSE  USAETL  1978 

39 

05 

02 

.22 

±0 

.20 

77 

18 

27 

.58 

±0 

.38 

LAYTON  1943 

39 

12 

45 

.52 

to 

.18 

77 

08 

32 

.89 

1+ 

o 

.38 

MAYNE  USAETL  1978 

39 

1 1 

13 

.90 

+0 

.20 

77 

07 

04 

.71 

to 

.38 

MILL  1958 

39 

20 

II 

.29 

±0 

.20 

77 

22 

49 

.22 

±0 

.38 

RIFFLE  USAETL  1978 

39 

08 

52 

.88 

±0 

.20 

77 

17 

29 

.15 

±0 

.38 

SHERWOOD  USAETL  1978 

39 

08 

54 

.81 

±0 

.20 

77 

00 

56 

.85 

±0 

.38 

SPENCER  1969 

39 

07 

14 

,06 

+0 

.20 

76 

59 

04 

.65 

±0 

.38 

TABOR  1969 

39 

15 

06 

.70 

to 

.20 

77 

08 

37 

.35 

±0 

.38 

TOWER  USAETL  1978 

39 

03 

40 

.08 

±0 

.20 

76 

55 

10 

.99 

±0 

.38 

VADER  USAETL  2978 

39 

10 

11 

.11 

±0 

.17 

77 

16 

36 

.92 

to 

.38 

WELFARE  1970 

39 

03 

15 

.95 

to 

.17 

76 

51 

23 

.40 

±0  .38 
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TABLE  13.  Deflections  of  the  Vertical 


STATION 


n  a 

n 


AERO  USAETL  1978 

+4 

".25 

+0  ".  1 7 

+2 

”.9 

±0 

".29 

BRINK  1 969  AZ.  MK. 

+  1 

.56 

+0 

.19 

+  1 

.8 

±0 

.29 

CEDAR  HEIGHTS  1969 

+2 

.49 

+0 

.20 

+2 

.4 

±0 

.29 

CHEVY  USAETL  1978 

-1 

.35 

±0 

.20 

-5 

.1 

to 

.29 

DALE  1943  RM  1 

+4 

.49 

+0 

.20 

+5 

.0 

±0 

.29 

EIRE  USAETL  1978 

+5 

.15 

+0 

.20 

+6 

.4 

±o 

.29 

E REEWAY  USAETL  1978 

+2 

.12 

.to 

.20 

+1 

.8 

+0 

.29 

HORSE  USAETL  1978 

+  i 

.63 

+0 

.20 

+3 

.0 

±0 

.29 

LAYTON  1943 

+  1 

.50 

±0 

18 

-0 

.9 

±0 

.29 

MAYNE  USAETL  1978 

+0 

.61 

±0 

.20 

-1 

.7 

+0 

.29 

MILL  1958 

+5 

.23 

+0 

.20 

+6 

.2 

±0 

.29 

RIFFLE  USAETL  1978 

+0 

.64 

±0 

.20 

+2 

.5 

±0 

.29 

SHERWOOD  USAETL  1978 

-1 

39 

±0 

.20 

-5 

.5 

±0 

.29 

SPENCER  1969 

-0 

.31 

+0 

.20 

-4 

.8 

±p 

.29 

TABOR  1969 

+1 

.42 

+o 

.20 

-0 

.5 

+0 

.29 

TOWER  USAETL  1978 

-1 

.28 

±0 

.20 

-5 

.2 

±0 

.29 

VADER  USAETL  1978 

+2 

.39 

IP 

.17 

+2 

.1 

+0 

.29 

WELFARE  1970 

-1 

.58 

+0 

.17 

-7 

.3 

±0 

.29 

Expected  Accuracy  Between  Two  Stations.  The  expected  accuracy 
between  stations  can  be  derived  by  applying  the  law  of  error  propagation.  The 
meridian  deflections  of  the  vertical  components  ||  and  £2  of  the  stations  1  and 
2  have  standard  errors  oj  and  o2  respectively. 

The  standard  error  of  the  difference 

•*  =  h  -  *. 
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(64) 


is  obtained  by  the  special  law  of  error  propagation: 


(65) 


and  £2  are  quite  independent  of  each  other,  the  covariance  o,  2  =  0.  Therefore, 


o 


a 


(66) 


From  the  values  shown  in  table  13,  the  following  is  obtained: 

(od)£  =  ±  0".20 y/T  =  ±  0  ".29  (67) 

In  a  similar  way,  the  standard  error  of  the  difference  of  the  deflections  of  the  vertical 
component  in  the  prime  vertical  is  found 

(od)n  =  ±0’’.41  (68) 


COMMENTS  ON  LATITUDE  DETERMINATION  •  The  latitude  deter¬ 
mination  by  the  N(iS  modified  Sterneok  method  involves  two  symmetrical  circum- 
meridian  zenith  distances  made  in  the  direct  and  reverse  theodolite  positions  east  and 
west  of  the  meridian.  The  index  error  is  eliminated  and  the  need  for  the  zenith  point 
correction  is  avoided.  An  additional  salient  feature  of  the  index  error  is  the  elimination 
of  the  personal  error  of  star  bisection. 


CONCLUSIONS 


1.  The  Sterncck  Method  of  determining  astronomic  latitude  entails  measure¬ 
ments  of  zenith  distances,  which  arc  subjected  to  circle  and  reading  errors, 
as  well  as  systematic  error  due  to  atmospheric  refraction.  These  errors, 
determined  by  a  relationship  between  zenith  distances  and  stars’  de¬ 
clinations.  are  excessive  for  the  needs  of  baseline  research. 

2.  Such  errors  are  difficult  to  control,  measure  or  eliminate  by  data  pro¬ 
cessing.  but  their  effects  can  be  minimized  by  a  special  least  square 
treatment  introduced  by  the  author. 

3.  A  new  observing  procedure,  associated  with  the  Wild  T  -  4  universal 
theodolite  and  observations  of  transit  times  of  pairs  of  stars  over  a  fixed 
vertical  plane,  significantly  reduces  the  number  of  unknowns  and  thereby 
provides  improved  precision  and  reliability. 


An  essential  condition  of  the  NOS  observational  procedure  is  that  the  in¬ 
strument  remains  clamped  in  the  meridian  and  the  star  is  observed  at  a  certain  distance 
from  the  middle  vertical  thread,  with  time  being  recorded.  TALC'OTT  level  vial 
sensitivity  values  are  highly  correlated  with  the  calibration  process  by  displacing  the 
level  division  intervals  with  respect  to  the  angular  change  of  the  vertical  circle. 


Absolute  zenith  distances  are  to  be  used  for  the  latitude  computation,  which 
is  hard  to  obtain  with  a  minimum  source  of  errors.  Anomalous  variations  in  the  re¬ 
fraction  affect  differently  the  north  and  south  slurs  and  the  error  in  level  calibration, 
which  increases  both  accidental  and  systematic  errors  of  observations. 


The  observer  probable  error  of  one  observation  is  0  '*.40,  This  small  probable 
error  is  a  criterion  of  the  skill  of  the  observer.  The  stations  latitude  mean  error  is 
less  than  0".::. 
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SYMBOLS 


X 

4> 

6 

a 

t 

Z 

P 

P 

e 

T 

A 

R 

i 


Longitude 
Latitude 
Declination 
Right  assension 
Hour  angle 
Zenith  distance 
Persona]  equation 

Level  sensitivity 
Index  error 
Interval  of  time 
Azimuth 
Refraction 
Inclination 
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